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Momentum widths of the primary fragments and observed final fragments have been investigated 
within the framework of an Antisymmetrized Molecular Dynamics transport model code (AMD-V) 
with a sequential decay afterburner (GEMINI). It is found that the secondary evaporation effects 
cause the values of a reduced momentum width, oo, derived from momentum widths of the final 
fragments to be significantly less than those appropriate to the primary fragment but close to those 
observed in many experiments. Therefore, a new interpretation for experiemental momentum widths 
of projectile-like fragments is presented. 
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Measurement of the momentum widths of projectile- 
like fragments (PLF) emitted in relatively peripheral nu- 
clear collisions have long been viewed as a means of ob- 
taining detailed information on the intrinsic Fermi mo- 
mentum distribution in the nucleus (l], [| ||, ||, ^, |t], 
|, | [tC], |ft], [TJ, f0| [Uj. This information is fundamen- 
tal in testing theoretical models of the nucleus and, in 
recent years, has taken on further importance in calcu- 
lations related to the production of secondary fragmen- 
tation beams now being employed to study both struc- 
ture and reactions far from stability fl5[| . Such studies, 
which explore a much broader range of neutron to proton 
asymmetry than previously accessible, are providing new 
insights into nuclear structure and nuclear astrophysics 
and offer the possibility of probing the nature of nuclc- 
onic matter in much greater detail fl6|| . 

In this Letter, we re-examine the validity of the 
Goldhaber Model |l|] within the framework of a non- 
equilibrium transport model, namely the AMD-V model 
of Ono et al JlTj] and study the effects of sequential decay 
on the observed momentum widths of final-state spec- 
tator fragments with the help of an evaporation code 
GEMINI @. We find that the Goldhaber Model works 
well in describing the momentum widths of the emerging 
primary fragments produced in the spectator fragmenta- 
tion. We find, further, that reduced momentum widths 
derived from the observed momentum widths of the final 
fragments are narrower than the reduced widths charac- 
terizing the primary fragment momentum distributions. 
This reflects the influence of the light particle evaporation 
from the primary fragments. Applying the Goldhaber 
Model without taking this evaporation contribution into 
account leads to an erroneous interpretation regarding 
the primary fragmentation step. 

In the pioneering work of Goldhaber Q, it was sug- 
gested that the momentum widths of PLFs are deter- 
mined by the intrinsic Fermi motion of the constituent 
nucleons which are removed from projectile during the 
break-up process. Assuming that (Aq — K) nucleons are 



suddenly removed from a nucleus which originally has Aq 
nucleons, a nucleus consisting of K nucleons will emerge. 
For this fragment of K nucleons, Goldhaber showed that 
the momentum width, a, could be expressed as : 

I k(a -k) 

a = a i Ao-1 (1) 

where oo, the reduced momentum width, is related to 
the intrinsic Fermi motion of a single nucleon. If the 
projectile nucleons have a mean square momentum in the 
projectile frame equal to (3/5)Pp, where Pp is the Fermi 
momentum, then that a momentum dispersion with <Tq 
— (1/5)P F is expected. Based upon electron scattering 
measurements of Fermi momenta |fl9f cto is expected to 
be~112-116MeV/c 0. However, experimental results 
favor values of a ~ 90 MeV/c §, | g §. 

In an attempt to explain this difference, Bertsch H 
treated the correlations between the momenta of indi- 
vidual nucleons localized in space and calculated a cor- 
rected reduced width for 40 Ar fragmentation which was 
17% smaller than Goldhaber's prediction. Murphy (t| 
considered the fact that the fragment is also a Fermi gas, 
ic., that there also exists a phase space constraint on the 
nucleon momenta in the projectile fragments which was 
not taken into account in the Goldhaber Model jj]], and 
predicted fragment momentum distributions which are 
narrower than those observed. In a purely kinematical 
semi-classical model Gan et al. B showed that the mo- 
mentum width of a fragment suddenly broken off from a 
Fermi distribution is sensitive to the single-particle distri- 
bution but the quantitative difference with the observed 
values still remains. Recently, Beshu et al. took into 
account the dependence of (Tq on the apparent tempera- 
ture of the PLFs, which are excited during the collision. 
In this case, the projectile fragmentation process is as- 
sumed to be a slower process and the width governed 
by the Fermi distribution of nucleon momenta in the ex- 
cited projectile as discussed by Bauer ||. In contrast to 
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most results for longitudinal momentum widths, Brady 
et al. observed broader transverse momentum widths 
than those predicted by the Goldhaber Model (l(], pi , 
These were explained by a collective motion: a bounce- 
off which imparts an additional transverse momentum to 
the spectators and results in wider momentum widths, 
especially for heavier projectiles. 

In recent years, several attempts have been made to 
provide a systematic phenomenological description of the 
available data ^, [l2], [lj| and the basic theoretical models 
have been extended S. In most of the theoretical work 
on the momentum widths of PLFs, analytical techniques 
have been employed to explore the problem. One paper, 
that of Gossiaux et al. [jlij , has employed a QMD model 
calculation to explore the properties of spectator matter 
in 600 MeV/nucleon collisions of Au with Au and con- 
cluded that the underlying physics is more complicated 
than that generally assumed in analytical approaches. 

The initial stage of the reaction studied in this paper, 
200 MeV/u 40 Ar + 27 A1, has been simulated using the 
AMD-V model AMD models have been very suc- 

cessfully used to study the static nature of light nuclei 
pjj . In this model, calculation of the minimum energy 
states can be carried out to define initial ground states 
of the colliding nuclei. 

In the AMD-V model the stochastic branching pro- 
cess of the wave packet diffusion is incorporated. The 
widths of the wave packets in each stochastic branch are 
kept constant and the dynamics of the widths of the wave 
packets calculated by the Vlasov equation is incorporated 
as a stochastic diffusion process of the centroids of the 
wave packets. This model was successfully applied to 
multifragmentation events in the 40 Ca + 40 Ca reaction 
at 35 MeV/u @ and 64 Zn + 58 Ni reactions Q at 35 - 
79 MeV/u. In the present calculations, the Gogny force 
p4j , which gave the best fit in these previous analyses, 
was used as an effective interaction. This Gogny force 
gives a momentum-dependent mean field and an incom- 
pressibility of 228 MeV for infinite nuclear matter. We 
note further that in our initialization of the projectile 
and target ground states in the AMD-V model, we get 
an initial width of the nucleon momentum distribution 
of - 105 MeV/c for 40 Ar projectile and - 101 MeV/c 
for the 27 Al target. These values are slightly less than 
expected from electron scattering measurements of Fermi 
momenta, ie. 112 - 116 MeV/c |lj|, but are still in the 
reasonable range. 

The calculation was started at t = 0, with a distance 
of 15 fm between the centers of the projectile and target 
in the beam direction. Each event was followed up to t 
= 300 fm/c. However, considering that we are initially 
interested in the properties of primary fragments, we an- 
alyze the AMD-V results at an earlier time, when the 
two dominant spectators just separate from each other 
and before significant evaporation can occur. 

Calculations were carried out in the peripheral collision 



zone, ie. b = 6 - 8 fm. 10000 events are calculated. Pri- 
mary fragment masses and their excitation energies were 
extracted at t s = 60 fm/c. Fragments were identified 
using a configuration space coalescence technique with a 
coalescence radius of 5 fm, but the size of the fragments 
depends only slightly on the coalescence radius. 

The primary fragments which are isolated in this way 
are excited. In order to simulate experimental data as 
closely as possible, a modified version of the GEMINI sta- 
tistical model code [|l8| has been used as an afterburner 
to follow the de-excitation of these excited fragments. In 
this modified version, discrete levels of the excited states 
of light fragments with Z < 14 are taken into account 
and the Hauser-Feshbach formalism is extended to treat 
the particle decay from a parent nucleus up to 50 MeV 
excitation energy. In the calculation, care was taken to 
follow the entire de-excitation cascade so that the distri- 
bution of initial parent nuclei (primary fragments from 
AMD-V) leading to each observed final fragment could 
be derived and the contributions of nucleon evaporation 
and mass loss to the momentum widths of final fragments 
could be determined. 

In order to study the primary fragments of spectator 
fragmentation, we first defined, in each event, the heavi- 
est fragment with a positive parallel velocity in the cen- 
ter of mass as a projectile-like fragment and the heaviest 
fragment with a negative parallel velocity in the center of 
mass as a target-like fragment (TLF). We then analyzed 
the three components of the momenta of the PLF and 
TLF fragments. Since the the collisions are peripheral, 
the primary PLF and TLF can be cleanly separated. 

These primary fragment momentum distributions can 
be described by a Gaussian shape characterized by a 
width a. The primary fragment mass numbers and the 
three momentum widths, ap x , ap y and ap z (Z is in the 
beam direction and X is in the direction of the impact 
parameter axis) are plotted in Fig. [l] as a function of 
the PLF and TLF fragment masses at b = 6-8 fm. The 
widths of the primary fragments are denoted by open 
circles and the dashed lines represent the results of fits 
to the Goldhaber Model expression (Eq.(l)). In order 
to obtain better fits, we allow both the initial mass Ao 
and the reduced width, <7q to be free parameters. The 
fit parameters are shown in upper-right corner. The ini- 
tial masses are close to the original masses of projectile 
or target and the average of the reduced widths of the 
three momentum components is ~ 105 MeV/c for PLFs 
and 110 MeV/c for TLFs. The values of <tq are very 
close to the initial nucleon Fermi momentum widths of 
the AMD-V ground states as mentioned before. Conse- 
quently, it appears that the Goldhaber model works well 
for the primary fragments. 

However, the calculation indicates that, the primary 
fragments are excited: the mean value of excitation en- 
ergy is about 1.7 MeV/u and 2.7 MeV/u for PLFs and 
TLFs, respectively. To evaluate the effect of secondary 
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FIG. 1: The three momentum widths crp z (upper), crp y (middle) 
and ap x (lower) as a function of the mass of PLF (left) and TLF 
(right) fragments. The open circles depict the momentum widths 
of the AMD-V primary fragments as a function of the mass of 
primary fragments and the dashed lines indicate the Goldhaber 
Model fits. The solid squares show the momentum widths of final 
fragments as a function of the mass of final fragment and the solid 
lines represent Goldhaber Model fits to those results. The open 
squares represent the momentum widths of final fragments as a 
function of the reconstructed average mass of the primary fragment. 
The solid circles depict the reconstructed momentum widths as a 
function of the reconstructed average mass of the primary fragment. 
See details in text. 



particle emission on the masses and momentum widths 
we used the AMD-V results at t s as input to the GEMINI 
calculation. As a result of the subsequent statistical de- 
excitation the final mass distribution becomes broader 
but the PLF and TLF components can still be identi- 
fied and the momentum widths for PLF and TLF final 
fragments can be determined. 

The solid squares in Fig. |l| represent the momentum 
widths of the final fragments, observed after the sec- 
ondary de-excitation. The results have also been fit with 
the Goldhaber Model expression (the solid line) and the 
extracted apparent primary mass Aq and width oo is 
shown on lower-left corner. The mass parameters remain 
close to the initial masses of the projectile and target 
but the cto values become significantly smaller, manifest- 
ing mean values ~ 83 MeV/c. Obviously, these values 
approach to the typical values of 90 MeV/c observed ex- 
perimentally (recall that the AMD-V ground state has a 
slightly lower Fermi momentum than those from exper- 
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FIG. 2: The primary fragment mass distributions for some selected 
PLF and TLF final fragments for 200 MeV/u 40 Ar + 27 A1 at 6 - 8 
fm. 



iments). It appears that the sequential decay plays an 
important role in causing the reduced momentum width 
of the final fragments to be narrower than that of the 
intrinsic nucleon momentum width. This leads us a rela- 
tively simple explanation for the experimentally reported 
momentum widths of PLF fragments. 

Why does the sequential decay make the reduced width 
decrease and not increase? In order to answer this ques- 
tion, we have tracked the secondary decay paths and 
obtained the distribution of parent nuclei (primary frag- 
ments) which lead to each of the observed final fragments. 
Fig. ^| shows some samples of the primary fragment dis- 
tributions which contribute to some selected final frag- 
ments. Clearly, the farther removed the final fragment 
from the primary parent nuclei, the broader the distribu- 
tion of primary fragments which contribute to that final 
fragment. 

If we plot the observed momentum width for the final 
fragments as a function of their average primary fragment 
masses, as represented by the open squares in Fig. [I], we 
immediately see that there is a large increase of momen- 
tum width in comparison with the original primary frag- 
ment width, denoted by open circles. It is, of course, 
natural that the sequential decay will make the momen- 
tum distribution wider. 

Knowing the parentage of a final fragment we can 
reconstruct the primary width distribution leading to 
that fragment. The reconstructed width for the av- 
erage parent fragments can be written as: (J recons = 
^/Yji-piO p r im(i) 2 j where pi is the fractional contribution 
of primary fragments to the final observed fragments. 
The results are indicated by the solid circles in Fig. [I]. 
As required by this procedure, the results fall on the 
curves defining the primary fragment momentum widths. 
The difference in widths represented by these points and 
those represented by the open squares reflect evaporation 
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FIG. 3: Mean reduced momentum widths for light particle emis- 
sion from PLF and TLF primary fragments for 200 MeV/u 40 Ar 
+ 27 A1 at 6 - 8 fm. 



effects. 

If we pursue the analogy to the Goldhaber Model, then 
from the decay calculations we can also define aQ Vap , a 
reduced width. This width characterizes the evaporation 
step and can readily be compared with the momentum 
width which results from the intrinsic nucleon momen- 
tum distribution from the primary step. We define the 
relationship between the width of the momentum distri- 
bution of final fragments of mass A b s and those of the 
primary fragments Ap r i m {i) which de-excite to produce 

A bs ■ 



a Aobs 



+ (A Prim (i) - A obs ){a* vap ) 2 )} (2) 



where pi is the fractional contribution of Ap r i m (i) to the 
final observed fragments (A b s ) and u^ vav is the contri- 
bution of particle evaporation to the momentum width. 
Fig. |3| presents the mean value of <jQ Vap over Px, Py and 
Pz directions as a function of the final fragment mass. We 
see that Ujf " ap is generally in the range of 30-70 MeV/c, 
significantly lower than Pp/^/E which characterize the 
primary break-up. 

For comparison, we also checked 40 Ar + 27 A1 at b = 
4-6 fm, the similar values and the same conclusion was 
drawn. In addition, we also checked our technique to ex- 
tract (To with the fixed Aq as is typically done in experi- 
mental analysis of initial projectile and target masses in 
the fit procedure with Goldhaber Model rather than the 
free parameter of Aq. The conclusions are not changed. 

In summary, the momentum widths of the primary 
fragments and final observed fragments produced in the 
reaction 200 MeV/nucleon 40 Ar + 27 A1 were investigated 
using an AMD-V calculation with a statistical model af- 
terburner. It is found that the Goldhaber Model works 
well for the primary fragments formed in spectator frag- 
mentation. The momentum widths of the primary frag- 
ments are basically related to the initial Fermi momenta 
in the projectile or target. However, since the primary 



fragments are excited, they are de-excited by light par- 
ticle evaporation. This secondary decay decreases the 
observed masses and increases the observed momentum 
widths of the primary fragments. As a consquence, it 
makes the reduced width parameter <jq, derived from a 
Goldhaber Model fit to the data, narrower by almost 20 
MeV/c than the initial Fermi momentum width, consis- 
tent with many experimental observations. Given that 
evaporation components can be experimentally deter- 
mined from particle- fragment correlation measurements, 
even in very complex reactions p5[ ], this result suggests 
that observed width distributions could be corrected for 
secondary decay and the Fermi momentum distribution 
of the primary fragment could be probed in relatively 
low intensity radioactive beam experiments. This would 
allow the extension of such measurements over a much 
larger range in N/Z ratio and provide even more strin- 
gent tests of our microscopic models of asymmetric nu- 
clear matter. 
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